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q | Abstract. In the present paper, we propose integrable semi-discrete and full-discrete 

D ' analogues of the short pulse (SP) equation. The key of the construction is the bilinear forms 

. and determinant structure of solutions of the SP equation. We also give the determinant 

formulas of /Y-soliton solutions of the semi-discrete and full-discrete analogues of the SP 
equations, from which the multi-loop and multi-breather solutions can be generated. In the 
continuous limit, the full-discrete SP equation converges to the semi-discrete SP equation, 
then to the continuous SP equation. Based on the semi-discrete SP equation, an integrable 
numerical scheme, i.e., a self-adaptive moving mesh scheme, is proposed and used for the 
numerical computation of the short pulse equation. 
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Q\\ 1. Introduction 
O 

Most recently, the short pulse (SP) equation 

X ; i , 

u xt = u+-(ur) xx (1-1) 

was derived as a model equation for the propagation of ultra-short optical pulses in nonlinear 
media HI 13- Here, u = u(x,t) represents the magnitude of the electric field, the subscripts 
t and x denote partial differentiation. Apart from the context of nonlinear optics, the 
SP equation has also been derived as an integrable differential equation associated with 
pseudospherical surfaces [3J. The SP equation has been shown to be completely integrable 
EH] [51 [6113. The loop soliton solutions as well as smooth soliton solutions of the SP equation 
were found in J8l [9l QU. The connection between the SP equation and the sine-Gordon 
equation through the hodograph transformation was clarified in HI], and then the soliton 
solutions including multi-loop and multi-breather ones were given by using the Hirota bilinear 
method. 

Integrable discretizations of soliton equations have received considerable attention 
recently [12l[T3l[TH[T5]|. In our recent work, the authors proposed an integrable semi-discrete 
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analogue of the Camas sa- Holm (CH) equation and apply it as a numerical scheme, i.e., a self- 
adaptive moving mesh scheme lfT6l[T7ll . The key of the discretization is an introduction of an 
nonuniform mesh, which plays a role of the hodograph transformation as in the continuous 
case. 

In the present paper, we attempt to construct integrable semi-discrete and full- 
discretizations of the SP equation by the same approach used in the CH equation. We also 
attempt to use the semi-discrete analogue of the SP equation as a self-adaptive moving mesh 
scheme to perform numerical simulations. 

The rest of the present paper is organized as follows. In Section 2, we review the bilinear 
equations and determinant solutions of the SP equation. In Section 3, we propose an integrable 
semi-discrete analogue of the SP equation, whose A'-soliton solutions are also constructed in 
terms of determinant form. By using the semi-discrete analogue of the SP equation as a self- 
adaptive moving mesh scheme, the numerical results for one- and two-loop solutions are also 
presented. In Section 4, the full-discrete analogues of the SP equation are proposed. The 
paper is concluded by Section 5. 

2. Bilinear equations and determinant solutions of the short pulse equation 

In this section, the results in |fTT| regarding the bilinear equations and the solutions of the SP 
equation will be briefly reviewed. 

First, by introducing the new dependent variable 



r 2 = l+u 2 x , 



(2.1) 



the SP equation is rewritten as 




(2.2) 



Introducing the hodograph transformation 

1 2 

dy = rdx + -u rdt , ds = 



dt , 



(2.3) 



i.e., 



d_ 

dt 



2 dy ds' 



i 2 a a 



_a_ 

dx 




a 



we obtain 




(2.4) 



where 



r 2 = 1 + r 2 ul . 



The equation (12.41) can also be cast into a form of 




(2.5) 
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Introducing new variables 
1 



u = h, (2.6) 



COS(|) 

Eq. (|2.5l ) leads to the sine-Gordon equation 

§ ys = sin((). (2.7) 
Moreover, as is shown in lfT8~ll2~2l . upon the dependent variable transformation 

^,,) = 2iln^4, 

the sine-Gordon equation (12.71) leads to the following bilinear equations 

1 
4 



FF ys -F y F s = ^{F 2 -F* 2 ), (2.8) 



F*F* S - F*F* = - (F* 2 —F 2 ) , (2.9) 

where F* is the complex conjugate of F. Henceforth, the solutions of the SP equation are 
obtained by F and F* through the dependent variable transformation 

„ (y , s)= ^ M = l( 2ito ^). (2 ,0) 

In what follows, we will show that the bilinear equations (I2.8I) -( |2.9I) are actually obtained 
as the 2-reduction of the two-dimensional Toda lattice (2DTL) equations: lfT9l 1201 12T1 1221 
1 2 

-DyD S Xn ■ X„ = X„ -T n +lT„_l , (2.11) 

i.e., 



d 2 X n d% n dl n 2 



"dYdS dY dS 
where D x is the Hirota D-operator which is defined as 



x„+iT„-i, (2.12) 



D "^=(sr!)" /WgWI -' 



Applying the 2-reduction % n -\ = a x n +i (a is a constant), we obtain 

3 2 T W dx n dx n _ 2 2 

"dYdS dYdS~ n n+u 1 ] 
where the gauge transformation x n — > QP-% n is used. Letting / = To and / = Ti, we have 

ff Y s-f Y fs = f 2 -f, (2.14) 

ffys-fyfs = f 2 -f 2 . (2.15) 
Under the independent variable transformation v = 2Y , s = 25, we obtain 

ffys-fyfs = \(f 2 -f 2 ), (2-16) 

ffys-fyfs^lif-f), (2.17) 
which are bilinear equations of the SP equation. 
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Next, we give the Casorati determinant (Af-soliton) solution of the SP equation. It is 
known that the Casorati determinant solution of the 2DTL equation is of the form ETll22ll : 



x„(y,5)= ^ +j - 1 \y : s) 



(2.18) 

l<i,j<N 



where (Y, S) satisfies linear dispersive relations 

For example, a particular choice of \]/^ (i 7 , 5) 

^ (Y, S) = c iAP y i¥+ * S+ ^ + a^e qiY+ k S ^ , (2.20) 

with and c ; .2 being constants, satisfies the linear dispersive relations and gives the N- 
soliton solutions. 

Applying the 2-reduction = — p, and the change of variables y = 2Y and s = 2S, we 
obtain the determinant solution of bilinear equations (12.161) and (12.171) : 

f(y,s) =% (y,s), f(y,s) =%i(y,s), 

(2.21) 



*n(y,s) = ^ n+J 1} (y,*) 



l<i,j<N 



where 

¥i W (y,,) = c i , l rie^*^ + cv(-iH) H e-^b'^ . (2.22) 

Since w is real and the dependent variable transformation w includes the imaginary 
number, we must consider the reality condition of u. Let us introduce a and [3 such that 
F* = af and F = (3/, where F and F* are complex conjugate of each other. Note that F and 
F* also satisfies the bilinear equations (12.161) and (12.171) because of 

■■sM-sM-s(»W-sH' »*> 

Thus a set of F and i 7 * gives solutions of the SP equation as well as a set of / and /. By 
choosing phase constants appropriately, the functions / and / can be made to be complex 
conjugate of each other to keep the reality and regularity of u. For example, the following 
choice 

= p n e jPiy+Wi S+ ^ 0i - in / 4 + (_ p .)" e -3^-2^+<- +i K/4 ? (2.24) 

guarantees the reality and regularity of the solution. 

Summarizing the above results, the determinant (/V-soliton) solution of the SP equation 
is given by 

x = y-2i(ln//),, t = s, 

f(y,s) =x (y,s), f(y,s) = %i(y,s) , 
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l<i,j<N 



where 



3. An integrable semi-discretization of the short pulse equation and numerical 
computations 

Based on the above fact, we construct the integrable spatial-discretization of the SP equation. 
Consider the following Casorati determinant: 



x„(M) 



l<i,j<N 
in) 

where \y\ satisfies the dispersion relations 

. (n) (n+1) 



-> (») 

W = -Ti 

Here A^ is the backward difference operator with the spacing constant a 

/(*)-/(*-!) 



>-i) 



(3.1) 

(3.2) 
(3.3) 



A k f(k) 



a 



Particularly, one can choose 

\fc\k,S) = C/ ,itf(l -a Pl )- k e^ S+ ^ + c l!2 q^l -aq^e^ , (3.4) 

which automatically satisfies the dispersion relations (13.21) and (13.31) . The above Casorati 
determinant satisfies the bilinear form of the semi-discrete 2DTL equation (the Backlund 
transformation of the bilinear equation of 2DTL equation) Il2~3"ll2~2l 

-D S - 1 ) X n (k + 1 ) • % n {k) + X n+1 (k + 1 )T n _i (*) = . 



(3.5) 



Applying 2-reduction 
and letting 



q -v = -Pi ■ 



N 



fk = to(k), fk = ti(k)= [Y[p 2 i U-i(k) 



d=l 



we obtain 



-Dsfk+l ■ fk - fk+lfk + fk+lfk = • 

a 

-Dsfk+l ■ fk — fk+lfk + fk+lfk = • 

a 



(3.6) 
(3.7) 
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where the gauge transformation % n — > (niLi^i) ^ 1S used. Note that / and / can be made 
complex conjugate of each other by choosing the phase constants properly. Under the change 
of independent variable s = 25, Eq. (|3.5l) implies the following two bilinear equations 

-D s f k+l -f k -f k+1 f k + f k+ J k = 0, (3.8) 
a 

-Dj k +\ ■ fk-fk+ifk + fk+ifk = 0, (3.9) 
a 

which can be readily shown to be equivalent to 

^ln^) _ 1 ) = fk±lA ) (3 . 10) 



\ a V fk J s J fk+lfk 

-(-(lnfy) -l)=M (3.11) 

\ a \ fk J s J fk+lfk 

Subtracting the above two equations, one obtains 

-(Ufy) -L^))=§±4-^. (3.12) 

a \ \ Jk J s V Jk J S J Jk+lJk Jk+lJk 

Introducing the dependent variable transformation § k ( s ) = 2iln (J^j^ 1 > one arrives at 

^ + i,-^ =s . n ^ + i+^V (3.13) 

which is nothing but an integrable semi-discretization of the sine-Gordon equation. Note that 
this is also known as the Backhand transformation of the sine-Gordon equation |[24ll25l . 
It is obvious that, from the semi-discrete sine-Gordon equation (13.131) . the equation 

( cos (^)) r _fcC, (3 ,4) 

is implied. By introducing the dependent variable transformation 

= if = 2 ""(awJ/ 8 >= ots (^H- <3 ' 15) 



it then follows 

"\ u k+i~ u k 



dhi- uh.-ui 



ds 4 

which is the first equation of a semi-discrete analogue of the SP equation 
From the facts 



(3.16) 



cos 2 ( *»\ + **) W ( *™ + *'\ = l , (3.17) 



2 J V 2 
sin = , (3.18) 



and 



2 J 2a 



1 §k ( 4>*+i+<hfc , 

— = — = cos , (3.19) 

r k a \ 2 
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it follows 



St {uk+\-u k ) 2 



+ ■ 



i.e.. 



4a 2 



2 (Uk+l—Uk)' 



(3.20) 



which becomes another equation of a semi-discrete analogue of the SP equation. 

Summarizing the above results, we obtained an integrable semi-discrete analogue of the 
SP equation and its solutions 



(wjt+i -u k ) 2 = 4(a 2 -8|), 

dd k 
ds 



2 _ 2 
U k+l ~ U k 



(3.21) 
(3.22) 



where the x-coordinate of the fc-th lattice point is given by X k = Xq + LrZo From the 
construction, the semi-discrete analogue of the SP equation has the following Casorati 
determinant solution: 

fk 



u k (s) 



d_ 

ds 



2iln 



fk 



a ( fk+lfk fk+lfk 



2 \ fk+lfk fk+lfk 



(3.23) 



k-l 
1=0 



fk(s)=%o(k,s), f k (s)=xi(k,s), 



where \|/- (k, s) satisfies 



\fi n, fcs) =p n l {\ -ap t )- k e*f ' ™ ■ + i 



x„(k,s) 



TS+^-o-iJt/4 



l<i,j<N 



'\+a Pl )- k e-^ s+ ^ +in,A . 



and the phase constants ±i7t/4 play a role of keeping the reality and regularity. 

Note that a 2 must be always greater than or equal to 8| because {uk+i —u k ) 2 > 0. This 
can be easily verified by 

'§k+i+§k 



\h\ 



a cos 



< \a\ 



(3.24) 



The mesh size of self-adaptive mesh |8jt| is always chosen as less than \a\. 

We can rewrite the semi-discrete SP equation in an alternative form which converges to 
the SP equation in the continuous limit S# — > 0. Multiplying Eq. (|3.22l) by 28^, we have 

^ = -5 t 

ds L 
Eliminating 8? using Eq. (l3.21l) . this leads to 



l k+l 



d(uk+i — Uk) 
ds 



§k(Uk+l+Uk) 



(3.25) 



(3.26) 
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Since 



d ( Uk+i~Uk\ _ 1 d(uk+i —Uk) Uk+i—Ukddk 
ds \ dk J $k ds 8? ds 



(3.27) 



it follows that 



d (uk+\—Uk\ Uk+\ +Uk fuk+i — Uk\ 

ds[—^) =Uk+l+Uk + ^^[—^) ' (3 - 28) 

by using Eqs. d3.26l) and (13.221) . Equation (13.281) gives another form of the semi-discrete SP 
equation. In the continuous limit a — > (8^ — > 0), we have 



u k+l —Ufr^du Uk+l + Uk 

8t dx ' 2 



u . 



dx_dx k y l d8j_ I* 1 , 2 1 2 

3* = 3/ + ^3* -> 3? - 7" 2 3t , 
os 4 



Consequently, Eq. (l3.28l) converges to 



3f - ^M 2 5x j = 2i/ + -MM 2 .. 



W + «(w x ) 2 -|--W 2 W xx , 



By the scaling transformation 2x — > jc, one arrives 

1 

2 l 

which turns out to be the SP equation 

«xf = U + -(U ) xx . 
O 

In a similar way employed in lTl6l[T7l . the semi-discrete analogue of the SP equation can 
be used as a novel numerical scheme, i.e., the so-called self-adaptive moving mesh method, 
to perform numerical computations for the SP equation. However, the first equation (|3.21l) 
has ambiguity for determining the sign even if the non-uniform mesh 8& is solved from 
the second equation (13.221) . To avoid this difficulty, we introduce an intermediate variable 
<j>£ = (4»jfc_|_i +<])fc)/2, and employ the following scheme, 

(uk+i-Uk) =2asin(<y, 

dh _ Uk+i +u k ^ 329 " > 
ds 2 

which can be derived from Eqs. d3.18l) and (13.151) . Equations (13.291) are equivalent to the 
integrable semi-discrete analogue of the SP equation, and the relation between the non- 
uniform mesh 8;t and <j>£ is 8^ = acos(^). Figures [T] and [2] are numerical results for one-loop 
and two-loop soliton solutions, respectively. The time stepsize is At — 0.01 and the number of 
grid points is N = 200. The detailed numerical results by using the integrable semi-discrete 
SP equation will be reported somewhere else. 
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Figure 1. Numerical solutions for one-loop soliton solution with (a) t = 0.0; (b) t = 10.0. The 
parameters of the initial condition are p\ =0.5 



4. Full-discretizations of the short pulse equation 

To construct a full-discrete analogue of the SP equation, we introduce one more discrete 
variable I which corresponds to the discrete time variable. 
It is known that the x-function 

x n (k,l) = ^ +J - l \kJ) , (4.1) 



l<i,j<N 



with 

^ n) (kJ)= Ci ^(l~a Pl )- k [ l-b- ) 



S+km + c ia q?(l-a qi )- k (l-b±-\ e^ 



satisfies bilinear equations [23 

'2 



-D s - 1 ) x n (k+ 1, /) ■ x n (k, I) + x„+i (k + 1 , l)x n -i(k, l)=0 : 



and 



Af-soliton solution, we obtain 

'2 



-D s - 1 )X n (k + 1, /) • X„(£, /) + X n+] (fc + 1 , OVt-1 0=0; 



and 



(2M), - 1)t»(M + 1) • t n +i(k,l) + x n {k, l)x n+ i{k,l + l)=0, 
where the gauge transformation x n — > (Y[f=i Pi) x n is used. Letting 

fk,i=x (k,l), fu=Xi(k,l), 



(4.2) 



(2bD s - l)x n (fc, / + 1) • x M+ i (*, /) + x n (k, l)x n+1 (k,l + l)=0. (4.3) 
Applying the 2-reduction x„_i = (njLi />; ) tn-t-l* i-e., adding constraints g ; = — p,- to the 



(4.4) 



(4.5) 
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Figure 2. Numerical solutions for the collision of two-loop soliton solution with (a) t = 0.0; 
(b) t = 6.0; (c) t = 8.0; (d) t = 10.0; (e) t = 15.0. The parameters of the initial condition are 
pi = 0.5,/j 2 = 1.0. 



the bilinear equations (14.41) and (14.51) imply the following four equations 

~ D s~ M fk+i,l ■ fk,i + fk+i,lfk,l =0, (4-6) 
-D s - 1 ) fk+u-fu+fk+ufu =0, (4.7) 



a 

(2bD s - \)f k . l+l ■f k ,l + fk,lfk,l+l =0, (4.8) 
(2bD s - l)/ w -fu+fufu+i =0, (4.9) 
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which are actually equivalent to 



z [x Jj±lA _ 1 + i±M&L = 0i ( 4.io) 



« V fkj J s fk+ljfkj 

^f {n fk+iA _ 1 + A±iA =0 , (4.11) 



a V fkj J s fk+ljfkj 
A fkj+i\ 1 , fk,lfk,l+l 



2b ln^ -l + ;'" '7 =0, (4.12) 

V JJU / v JkJ+lJkJ 



2b 



L!m±i) _ l + fk,ifu + i =Q (413) 

V A,/ / s fk,l+lfk,l 



Note that / and / can be made complex conjugate of each other by choosing the phase 
constants properly. By introducing 

u k ,= (nin^-) , (4.14) 



fkj 

and 

X k j = ka- (In fuf k: i) s , (4.15) 
where X k i is the x-coordinate of the fc-th lattice point at time /, we find the following relations 

. ( fk+l,lfk,l fk+l,lfk,l\ iA 

Uk+i,i - u k ,i = ia -= =— , (4.16) 

\Jk+l,lJk,l Jk+lJJkJ/ 

1 f fkjfkj+l fk,lfk,l+\\ ,a m\ 

Uk,i+i+u k ,i = r [- = -= — , (4.17) 

v \Jk,l+lJk,l Jk,l+\Jk,lJ 

Xk+ll - XkJ = a -(pl4 L + § ±1 4 L ), (4.18) 

2 \Jk+l,lJk,l Jk+lJJkJ/ 

X k j + i-X u = -- + — - + ^ — ■ (4.19) 

o lb \jkj+ljkj Jk,l+lJk,U 

It is straightforward to derive 

{u k+hl -u ktl ) 2 = 4{a 1 -§l l ) 1 (4.20) 
from Eqs. d4.16ft and (14. 1 8b and 

If P 2N 



(»*./-!+»*./)" = 4 ( jp-\X k j +1 -X ki i + -J ). (4.21) 

from Eqs. d4.17l) and (14. 19b . where = X k+U -X k j. Equations (14.201) and (14.211) give a 
full-discrete analogue of the SP equation. 

Let us consider another full-discrete analogue of the SP equation. From Eqs. (|4. 161) - 
Eqs. d4.19l) . we obtain 

fk+ljfkj 1 f Y y j Uk+1,1 ~ Uk,l \ 

7 7-= -\ X k+l,l-Xk,l-l = , (4-22) 

Jk+LUkj a\ 1 J 

fk±Ufk£ l ( v y i : "*+M ~ u Kl \ rA0 ^ 

J — J - = -[X k+u -X k , l + 1 —J, (4.23) 
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b(x kJ+I -X u+ 1 T - i UtJ+ \ +UtJ ), (4.24) 



fk,l+lfk,l V ^ 

/m/m+i -^ Zfem _ Zu+ l +i M W + ^/ ) (425) 



fk,i+ifk,i \ b 2 

From the relations (I4.22I) - (I4.25I) . we have 

x k+l,l+\ ~ x k,l+l -1 9 



(4.26) 



v v ; u k+i,l— u k,l 
^k+lj—^kj— 1 2 

_ x k+l,l+l-^k+l,l + b~ 1 2 

V V I 1 I \ u k.l+\+ u k,l 

A k,l+l ~ A k,l + b +l ^^l 

Equating the real part and imaginary part respectively, we have 

tv v \fv v , 1\ , u k+U+l- u k,l+lUkJ+l+Uk,l 
[X k +l,l+l ~ X k ,l+l) I *k,l+l -AJt,/ + - I H 

/ y y , !\ ,y y \ U k+l,l+l+Uk+l,l U k+U -U k J 

= x k+l,l+l ~ x k+l,l + t K-X-k+1,1 - x k,l) ~ ~ j (4.Z/J 



= ^-Xjfc+l^+l — + ( u k+l,l ~ u k,l) + -Xk,l)( u k+l,l+l + u k+l,l) > (4.28) 

which can be rearranged into the following simpler form: 

(X k +l,l+l — Xk+l,l —Xk,l+1 +%k,l) (r —Xk+l,l+Xk,l+l 
_ _ U k +l,l+l + Uk+l,l — Uu+1 — Uk£ Uk+l,l + Uk,l+1 

- 2 2 ' ( y) 

— Uk+i,i — uu+i +uk,i) (- +X] c+ ij + \ —X k+ ij+X k j + \ —Xkj 

= (X k +i,i+i +X~k+\,i —X^i+i — -Xjt,/) +«jfc+i,z + Kjfc,j+i + • (4.30) 

Equations (14.291 ) and (14.301) constitute another form of integrable full-discretization of the SP 
equation. Taking the continuous limit b — > in time, we obtain 

(^jfc+l -X k ) s = -^{uk+l-Uk){uk+i + u k ) , (4.31) 

and 

{u k+ i-u k ) s = (X k+1 -X k )(u k+l +u k ) , (4.32) 

which are nothing but the semi-discrete analogue of the SP equation (|3.21l) and (13.221) . Here 
we used J± ^- J — » d s F as b — > 0. 

From the construction of the full-discrete analogue of the SP equation, the determinant 
solution of the full-discrete SP equation is 

. ( 8k,l 8k,l\ 3 / f k j\ ,. 
Xtj = * fl " I (ff + ff) = * fl " ^('"A'A,') . (434) 
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gk,i = Po(k,l), gu = Pi(M) 



x n (kj) 



where \\f^ (k, I) satisfies 
Mff ) {k,l)=p n i {l-ap i )- k 



^\k,l) ^ +1 \k,l) 
Vt 1] (k,l) V { 2 n+l \k,l) 



^(M) ^ +1 \k,l) 



v { r N - i \k,i) 
vr N - i \k,i) 



v { r N - i \k,i) 
v { r N - i] (k,i) 



v { r N - 1] (k,i) 



1-.1 
Pi 



1 



rS+nio+m/4 



+ (- Pl ) n (l+a Pl )- k ^l+b-j e ^ 

and the phase constants ±in / 4 play a role of keeping the reality and regularity, s is an auxiliary 
parameter. Note that p™ can be expressed as p™ = 2d s X n (k, I) because the auxiliary parameter 
s works on elements of the above determinant by 2d s \\f\ n \kJ) = (k, I). In the lattice 

KdV and lattice Boussinesq equations, one of x-functions is also expressed by the derivative of 
another x-function with respect to an auxiliary parameter Il26ll27l . This is a common property 
of discrete soliton equations which are directly connected to the Backlund transformations of 
continuous soliton equations. 

Let us consider Eqs. (14.20l) and (|4.21l) again. Rewriting Eqs. (|4.20l) and (14.211) . we have 

2 



Uk+lJ—Uu 



u k,l+l + UU 



+ d 2 =a 2 



+ X, 



v k,l+l 



X kJ + T) = £2 ■ 



(4.35) 
(4.36) 



These equations actually give conserved quantities because a 2 and l/b 2 are constants. 
Introducing 



Iu = 



Jk,I 



Uk+l,l — uil 



UkJ+\ +Uk,l 



+ X, 



^kj+l 



Xu + - 



1 



Eqs. (|4.35l) and (14.361) imply the following conserved quantities 
? 1 

b 2 ' 



hi = a , Jkj 



(4.37) 
(4.38) 

(4.39) 
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for arbitrary integer values of k and /. Hence, we have 

/*,/+!-/*,/ =0, J k +i,,~J k j = 0. (4.40) 
A substitution of the corresponding conserved quantities leads to 

Uk+l,l+l + Uk+i,i — Uk.i+\ — Uk,i \ / — ~ u k,i+i + H,l 



= —(Xk+i,l+i +Xk+i,i — -Xjfc,z+i — X^/) (Xjt+i^+i — Xjfe+i,/ — Xk,i+i +%k,i) j (4.41) 

W/t+l,/+l + W /U+1 + M Jt,A / M fc+l.Z+l + u k+l,l ~ u k,l+l — u k,l 



— ( Xfc+l.Z+l — -Xfc+i.z +Xu+l — -Xjfc,/ + 7 ) / + i — Xk+ij —Xkj+i +Xh{) . (4.42) 



It can be readily shown that the difference of Eq. (|4.42|) and Eq. (|4.41|) gives Eq. (|4.29l ), whereas, 
the quotient is nothing but Eq. (14.30l) . In summary, Eqs. (l4.35l) and (14.361) . which imply 
conserved quantities, can also be derived from the full-discrete analogue of the SP equation 
(14391) and (14301) . 



5. Conclusions 



In the present paper, we proposed integrable semi-discrete and full-discrete analogues of the 
short pulse equation. The iV-soliton solutions of both the continuous and discrete SP equations 
were formulated in the form of Casorati determinants, which include multi-loop soliton and 
multi-breather solutions. Based on the semi-discrete SP equation, a self-adaptive moving 
mesh method is proposed and used for the numerical solutions of the SP equation. The 
examples of one- and two-loop soliton solutions shows the potential of this novel method 
for the numerical study of the short pulse equation. 
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